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Abstract 

We investigate the cosmology of a class of model with noncanonical scalar field and mat¬ 
ter both in FRW closed and open background. Writing the Einstein Equations in terms 
of dimensionless dynamical variables suitable for studying bouncing solutions, a fixed point 
analysis is carried out. Cosmological solutions satisfying the stability and bouncing condi¬ 
tions are obtained. 


1 Introduction 

There are two scenarios exist in literature, namely, inflation and bouncing model which address 
the shortcomings of the standard model of cosmology. Though inflation solves most of the 
problems (horizon, flatness and entropy) of the standard model of cosmology, the issue with 
the initial singularity is not resolved under its domain. It is the alternate scenario, nonsingular 
bouncing model, that eradicates the singularity by constructing a universe which begins with a 
contracting phase and then bounces back to an expanding phase through a non zero minimum 
in the scale factor. Nonsingular bouncing models can be categorised into two types, matter 
bounce model [T] and and Ekpyrotik models [2]. Eor a review on these models refer to 
and [1]. 

In this paper we consider a noncanonical scalar field with a general function of kinetic term 
F{X), where X = These theories are originally motivated to provide a large tensor 

to scalar perturbation in inflationary settings [SHZj. Dark energy with a general kinetic term 
F{X) is modeled first in ref. [8]. For other variants of models of dark energy in this context 
refer to [9]. Other works related to unifying dark matter, dark energy and/or inflation for 
noncanonical scalar held models are studied in [101413) . In order to study the phase space in this 
model, we write the hrst order equations of motion in terms of dimensionless dynamical variables 
|14j . The motivation to use noncanonical scalar held as matter is to construct nonsingular 
bouncing models. The phase space analysis of a cosmological model with scalar held Lagrangian 
F{X) — E(</) and matter for an FRW background is given in ref. [15j . The condition for 
nonsingular bounce is also discussed in ref. m- In order to explore the behaviour of curvature 
parameter near bounce in an nonsingular bouncing model we do a phase space analysis in an 
FRW closed and open universe. This can be easily extended to other nonsingular bouncing 
models. 

We study the cosmology of an anisotropic universe with a matter Lagrangian of the form 
F{X) — V{4i) and matter. In section [2] we write the Einstein’s equation in terms of dynamical 
variables suitable for the analysis of a bouncing scenario in FRW closed and open universe. Then 
we discuss hxed points and their stability in section [S) Conditions for existence of nonsingular 
bouncing solution is derived in terms of dynamical variables in section H) We summarise our 
results in section E) 
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2 Einstein Equations in Bianchi I background 

The action for our model is given by 

S = J d^x^i^R + F{X) - V{cl>) + Lm] (1) 

where is the lagrangian of the matter field. 

To see the behaviour of curvature parameter of the spacetime in a nonsingular bouncing 
scenario we work with a FRW closed and open universe. The line element of the same is given 
by: 


ds^ = —dt^ + a^{t )\———^ + r‘^siv?‘6d<jP‘]. (2) 

1 — kr‘^ 

where k = +1 denotes closed and k = —1 denotes an open universe respectively. 

The Hubble parameter H is defined as 


In terms of Hubble parameter, the 


H = 


1 da 
a dt 


Einstein equations take the following form 


dH 

dt 


- ^(/3 + 3p), 

p k 

3 " 


where p = p^ + Pm and p = p<j,+ Pm- 

Here the energy density p^ and pressure of the scalar field is found to be 


( 3 ) 


( 4 ) 


P(j> = 2XFx — F + V, 

p^ = FiX)-Vi^), (5) 


and Pm and pm are the energy density and pressure due to the term Lm- 
Substituting Eq. ([S]) in first and third line of Eq. (jH) , we get 

^ - \{2XFx -F + V + Pm + ‘i{F-V) + ^Pm) 

dc 0 


( 6 ) 


2 ‘2XFx - F V k Pm 

-W = -^- \- — - K -\ - 


3 ' 3 ' 3 

Here we further define few more variables which are useful for defining dimensionless dy¬ 
namical variables. They are 


Pk = 2XFx — F , 

F 

~ 2XFx - F ’ 

1 dlogV 

V3\pk\ dt 


( 8 ) 


where pk is the kinetic part of the energy density p^p, Wk is the ratio of kinetic part of the 
pressure p^p to the pk and a is the auxiliary variable which depends on the variation of potential 
with time. 
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Neglecting the interaction between scalar field and matter, continuity equation for in 
terms of dimensionless time variable N [dN = Hdt), is 

^ {2XFx -F + V) + 6XFx = 0. (9) 

Now we define a set of dimensionless dynamical variables which is suitable for nonsingular 
bounce models. Relevance of these variables that they remain finite during the entire evolution 
across bounce. The dynamical variables are 


X = 


> y = , z = ^^^=sign{z ), Qrn = 

v\Pk\ V\Pk\ VlPk \Pk\ 


( 10 ) 


Here sign{z) =sign{k) denotes FRW closed universe for +1 and open for —1. Using Eas. lfTOl) . 
Q and Q and parameters defined in Eq.([5]), the evolution equations of x, y and z are written 
as, 


di 

dN 

dy 

dN 

dz 

dN 

dflrn. 

diV 


= “2 {Wk - Wm)sign{pk) + (1+ Wm)ix‘^ - y\y\) ++ ^Wm) 

3 

+ ^x[{wk + l)x - ay\y\sign{pk)], 

3 

= ^y[-(^ + {wk + f)x-cry\y\sign{pk)], 

3 

= -H+-z{x{l + Wk)-y\y\sign{pk)), 

= -3(1 + Wm)x^m - [3ay\y\sign{pk) - 3x(l + Wk)] , 


where dN 



dt and the constraint equation relating dynamical variables is 
x'^ - y\y\ + = 1 X sign{pk). 


( 11 ) 


( 12 ) 


The equation for parameter a becomes [K] 


dcr 

diV 


-3^2 (T - 1) + 


3(7 (2S {wk + I) + Wk - 1) 
2 (2c7 + 1) {wk + 1) 


{wk + 1) X - af 


(13) 


where - = and E = 

Eor our model we have taken power law form for F{X) = FqX^, where Fq is a constant. 
Eor this form of F{X), Wk = 2 i^ H = r/ — 1. 

Potential V{(j)) is taken as V{(j)) = where Vq and c are constants with positive values. 

Eor this choice of V{(f)), T becomes unity. 

In the next section, we do a fixed point analysis of dynamical equations for x, y, z and a. 
The evolution of flm is determined from the constraint Ea. (|12l) . 


3 Fixed Point Analysis 

In this section, we do a fixed point analysis of our system of dynamical equation in order to 
extract the qualitative information about the nature of solution. Eixed points are calculated by 
taking the first derivative of the dynamical variables to be zero. The stability of a fixed point 
is determined from the behaviour of a small perturbation around that fixed point. 
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We get the set of fixed points Xc, Vc, Zc and Uc by solving the following set of equations 
simultaneously (where the subscript c denotes fixed points). Now, if we define the slopes of the 
dynamical variables x, y, z and a as f{x, y, z, a), g{x, y, z, a), h{x, y, z, a) and i{x, y, z, a). The 
set of equations we need to solve to obtain the fixed point is 


f{x,y,z,a) 

g{x,y,z,a) 

h{x,y,z,a) 

i{x,y,z,a) 


dx 

-^ = 0, 

dN 

dy 

i = o. 

dN 

da 

-^ = 0, 
dN 


(14) 


where, 


f{x,y,z,a) 


g{x,y,z,a) 

h{x,y,z,a) 

i{x,y,z,a) 


3 ,21 I 

= “2 - Wm){signpk) + (1 + - y\y\) + (1 + 

3 

+ 2 ^[{wk + l)x - ay\y\sign{pk)] , 

3 

= ^yl-^ + + ^)x - ay\y\sign{pk)], 

= -3zx + 3zx{l + Wk) - 3zy\y\sign{pk), 

_ 3 [2r.{wk + 1) + {wk — 1 )] 

2 2{2a + l){wk + l) 


-[(rcfc + l)x - cry^]. 


(15) 


The corresponding fixed point for 0,^ can be found using the constraint Ea. (ll2p . 

The stability of the fixed points can be examined from the evolution of pertubations around 
fixed points. Now, if (xc, yc, Zc, ac) is a fixed point and 6x = x — Xc, Sy = y — ijc, 6z = z — Zc and 
6a = a — ac he the respective perturbation around it, then the evolution of the perturbation is 
determined by 


6x = X = f{xc + 6x,yc + 6y,Zc + 6z,a + 6a), 

Sy = y = g{xc +6x,yc +6y,Zc +6z,a + 6a), 

6z = z = h{xc + 6x,yc + 6y^Zc + 6z,a + 6a), 

6a = a = h{xc +6x,yc +6y,Zc +6z,a + 6a), 


The evolution equations,upto first order, for these pertubations are 


where the matrix is 


^ 6x 


(5x \ 

6y 

= A 

6y 

6z 


6z 

\6a ) 


\6a ) 


A = 


/ ^ 

dx 

dx 
dh 
dx 
di 
dx 


at 
\ dx 


dv 

dg^ 

dy 

dh 

dy 

m 

dy 


dz 

dg 

dz 

dh 

dz 

di 

dz 


da 

da 

dh 

da 

Mi 

da 


\ 

/ 


(16) 


(17) 


(18) 
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Fixed Points {xc,yc, Zc, cTc) 

Stability Conditions 

( 0 , 0 , 0 , 0 ) for Wk = Wm 
( 1 , 0 , 0 , 0 ) 

(- 1 , 0 , 0 , 0 ) 

(0,0, ) 0) with Wm > Wk and Wm > 3 

^ Wm> 1 

Can’t decide 

Stable (see Fig.pP) 
Stable (see Fig. [2]) 
Can’t decide 

Can’t decide 


Table 1: Stablity Analysis of fixed points for closed universe with sign{pk) = +ve 


is the Jacobian matrix and is evaluated at the fixed point (xc, Vc, <7c) and hence each entry 
of A is a number. The solution of the system of equations can be found by diagonalizing the 
matrix A. A non trivial solution exists only when the determinant |A — AI| is zero. Thus, 
solving this equation in A we would get all the eigen values of the system corresponding to each 
fixed points. 

We have two cases: one with positive kinetic term, sign{pk) = +ve and other one with 
negative kinetic term, sign{pk) = —ve.. 


3.1 Closed Universe 

3.1.1 Case I, with signpk = +ve 


In this case we study the hxed points for all possible values of parameters in an FRW closed 
universe. The hxed point (0,0,0,0) is obtained for = Wm signifying all the dynamical 
variables x, y, z and u, going to zero at late times. It is a nonhyperbolic hxed point as the eigen 
value of A for this is (0,0,0,0). It’s stability cannot be decided from our hrst order analysis of 
perturbations. From now onwards eigenvalues would mean eigenvalues of matrix A. 

The second hxed point (1,0,0,0) denotes a late time kinetic dominated universe with other 
dynamical variables y, z and cr becoming zero. In this case eigenvalues are ^ _i _j_ 

|(1 + Wk), |(—1 + Wk + 3{wk — Wm))- This is a stable hxed point for the region of 

parameter space shown in the Fig. [T]. 

The next stable hxed point in this subsection is (—1,0,0,0) with eigenvalue (1 + |(—1 — 


Wk), |(-l-t(;fc). 




{-1-Wk) — ^{1 + Wk) + 3{1 + Wm)) shows again 


2(1+I0fc) ’2V 2' 

a late time kinetic dominated phase but with a negative value of Hubble parameter H signifying 
a contracting universe. This hxed point is found to be stable for the region of parameter space 
shown in Fig. [2]. The point ( — 1,0,0,0) may not be important from bouncing point of view, 
as we need the universe to transit to an expanding phase to be discussed in section HI 

The remaining hxed points, in this section, being (0,0,0) for Wm > Wk and 


Wm > —V, ( 0 , ^ , 0 , 0 ) with Wk > Wm and Wm > — Ij with eigen values 

( 0 , 0 , 


+ ^wl - Wm- ^WkWm, \J\\J ^k + 2wl - Wm - SWkWm), ( 0 , 0 , -1 

'^\J'^k + w"^ — Wm — WkWm) are also nonhyperbolic points. The stability of such hxed points 
goes beyond the linear stabilty analysis. All the hxed points and their stability conditions are 
noted in tabl^H 


3.1.2 Case II, sign{pk) = —ve 


In this section, we state the results of stability anlysis of our dynamical variables for the negative 
sign of kinetic energy density. The hxed points are found to be (0,0,0,0),(0,0, —\/0) 

and ( 0 , ^ 0 , 0 ) with eigen values ( 0 , 0 , 0 , 0 ), ( 0 , 0 , -\f^\J-Wk - ^wl + Wm + SwkWm, 
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Figure 1: Allowed region of parameter space for the fixed point (1,0, 0,0) in closed universe 
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Figure 2: Allowed region of parameter space for the fixed point (—1,0,0,0) in closed universe 
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Fixed Points {xc, yc, Zc, Oc) 

Stability Conditions 

(0, 0,0, 0) for Wk = Wm 

(».o. 

Can’t decide 

Can’t decide 

Can’t decide 


Table 2: Stablity Analysis of fixed points for closed universe with sign{pk) = —ve 






Figure 3: Evolution of the dynamical variables x {top left), y {top right), z {bottom left) and 
a {bottom right) for the fixed point (xc, <7c) = (1)0,0,0) with the values of parameters 
sign{z) = +ve, Wk = —2.0, Wm = 1/3 and sign{pk) = +ve for different initial conditions 


, I-in, - 37;;2 -U -C i and iO 0 _o o .pstipc 

tively. All these fixed points are nonhyperbolic and tabulated in tabl^SJ 

3.2 Open Universe 

3.2.1 Case I, with signpk = +ve 

In this case we study the fixed points for all possible values of parameters in an FRW open 
universe. The fixed point (0,0,0,0) is obtained for Wk = Wm signifying all the dynamical 
variables x, y, z and a, going to zero at late times. It is a nonhyperbolic fixed point as the eigen 
value of A for this is (0,0,0,0). It’s stability cannot be decided from our first order analysis of 
perturbations. 

The second fixed point (1,0,0,0) denotes a late time kinetic dominated universe with other 
dynamical variables y, z and a becoming zero. In this case eigenvalues are ^ _x -g 

|(1 + Wk), |(—1 + Wk + ^ 3{wk — Wm))- This is a stable fixed point for the region of 
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Fixed Points (xc, Vc, Zc, (Tc) 

Stability Conditions 

(0,0, 0,0) for Wk = Wm 
(1,0, 0,0) 

(-1,0,0,0) 

(0, 0, ^ ^ijXwS ’ ^ Wm> 3 

(0, , 0,0) with Wk > Wm and Wm> 1 

Can’t decide 

Stable (see Fig. [3]) 
Stable (see Fig. [5]) 

Can’t decide 

Can’t decide 


Table 3: Stablity Analysis of fixed points for open universe with sign{pk) = +ve 



Figure 4: Allowed region of parameter space for the fixed point (1,0,0,0) in open universe 


parameter space shown in the FigUJ 

The next stable fixed point in this subsection is (—1,0,0,0) with eigenvalue (1 + |(—1 — 
3(-l-«)fc)(-l+«jfc+ ) 

Wk), K-l-tCfc),-2 (1+I»fc) - - -’ ^{-1 + Wk)-^{1 + Wk) + 3{1 + Wm)) shows again 

a late time kinetic dominated phase but with a negative value of Hubble parameter H signifying 
a contracting universe. This fixed point is found to be stable for the region of parameter space 
shown in Fig.|5]. The point (—1,0,0,0) may not be important for bouncing point of view, as 
again, we need the universe to transit to an expanding phase to be discussed in section HI 


The remaining two fixed points being (0,0, — 


\/1+3'ih77 


, 0) with Wk > Wm and Wm > — 


and (0, ^ with Wk > Wm and Wm > — 1 with eigen values 

(0,0,+ 2>wl - Wm - 3WkWm, \J\\JWk + ^wl - Wm - 3WkWm) and 

(0,0, --^\Jwk + w‘l-Wm- WkWra, -^\J Wk + w^ - Wm - WkWm) are also nonhyperbolic points. 
The stability of such fixed points goes beyond the linear stabilty analysis. All the fixed points 
and their stability are noted in table El 


3.2.2 Case II, sign{pk) = —ve 

In this section, we state the results of stability anlysis of our dynamical variables for the negative 
sign of kinetic energy density. The fixed points are found to be (0,0,0,0),(0,0, — 0) 
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Figure 5: Allowed region of parameter space for the fixed point (—1,0,0,0) in open universe 


Fixed Points {xcUci ZcCTc) 

Stability Conditions 

( 0 , 0 , 0 , 0 ) for Wk = Wm 

( 0 , 0 , ^ , 0 ) with Wm > Wk and Wm > 3 

(0, ^J Wm> 1 

Can’t decide 

Can’t decide 

Can’t decide 


Table 4: Stablity Analysis of fixed points for open universe with sign{pk) = —ve 


with Wm > Wk and Wm > “ji and (0, y ^(^^^,0, 0) with Wm > and Wm > — 1 with eigen 
values ( 0 , 0 , 0 , 0 ), ( 0 , 0 ,- 3rc| -|- Wm + SwkWm, \I\\J-Wk - 3rc| + Wm + SwkWm) and 

/n n -Wk-uii+uim+WkW'm N 1 All J-I, c J ■ J- 

(0,0,- ^^^ - ) respectively. All these hxed points are non- 

hyperbolic and tabulated in table. SI 


4 Bouncing Scenario 


Now we obtain the conditions for non singular bounce to occur and also show the evolution 
of dynamical variables numerically. A nonsingular bounce is attained whenever the universe 
passes from a contracting phase to an expanding phase through a minimum value of the avearge 
scale factor a(t), but not zero. Mathematically, it satisfies 


_ 1 / da{t)\ 


= 0 , 


ab{t) \ dt Jb 

where subscript b denotes value of the variable at the bounce, and 

/ d^o(t) 




> 0 


(19) 


( 20 ) 
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Figure 6: Evolution of the dynamical variables x {top left), y {top right), z {bottom left) and 
a {bottom right) for the fixed point {xc,yc,Zc,o'c) = (1,0,0,0) with the values of parameters 
sign{z) = —ve, = —2.0, Wm = 1/3 and sign{pk) = +ve for different initial conditions 
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for minimum to occur. This implies 

mr(S)AD> 

Now, writing the above conditions in terms of dynamical variables for bouncing, we get 
Xb = 0 and > 0 which translate to the following equation 



dx \ 3 


{Wk - Wm){signph) + (1 + Wm){-y\y\) + 


> 0 . 


.divA 

This implies 

(y\y\{l + Wm) - 5|^|(1 + > 1 x sign{pk){wk - Wm)- 

At the bounce we then obtain the constraint equation among dynamical variable as 


( 22 ) 

(23) 


- y\y\ + zjzj - = -y\y\ + z\z\ - = 1 x sign{pk). (24) 

Now, for different negative initial conditions of x (contracting phase). Fig. [3] {top left) 
and [6] {top /e/t)for closed and open universe respectively, show its transition to positive values 
(expanding phase) crossing zero (bounce). The bouncing is guarenteed by the positivity of the 
slope of X as shown in Fig. [7] (left plot for closed and right plot for open). Thus, top left of Fig. 
[3] and Fig. [6] and Fig. [7], together, do indeed represent stable bouncing scenario in FRW 
closed and open universe. This is obtained by setting the values of equation of state parameters 
Wk = —2 (r/ = 1/4), Wm = 1/3 and sign{pk) = +ve and sign{y) = +ve. The evolution of other 
dynamical variables can be seen in Fig. [3] and [3], which show their asymptotic evolution to 
the respective fixed points for the same choice of parameters. 

It can be seen that the fixed point {xc, yc, Zc, <7c) = (1,0,0,0) does give rise to a stable bounc¬ 
ing universe as it satisfies Eas. (f23]l and [24)1 for open {sign{z)=-ve) and closed {sign{z)=+ve) 
universe. From this analysis, we conclude that finally after the bounce our universe at late 
times is driven by kinetic energy density in both the cases. The other fixed point (—1,0,0,0), 
though stable, can not give rise to a bouncing scenario as it ends up with a negative value of 
Hubble parameter, H, signifying a late time contracting phase. 

Also, we show the behaviour of curvature parameter, z, in this nonsingular bouncing set up. 
The curvature parameter increases initially in the contracting phase reaching an extremum at 
the bounce and then decreases to zero in the expanding phase as shwon in Fig. [8] for both open 
and closed universe. Thus, the curvature parameter remains finite at the bounce as expected in 
a nonsingular bouncing scenario and at late time universe becomes flat irrespective of whether 
we start initially with closed or open.This may be useful for building realistic models. 

The comparision between bouncing solutions for open and closed is done in Fig. [9]. It 
has been found that bouncing occurs earlier in the case of open than in the closed universe as 
shown in left hand side of Fig. [9] for the same set of initial conditions and parameters. Also, 
it is noted that, though, the solutions differ appreciably near the bounce, they approach to 
the same value at late time owing to zero value of the curvature parameter. The nonsingular 
bounce happens only for negative values of with our choice of parameters as shown in Fig. 
[9] {right) for both open and closed universe. 

Finally we show the effect of different values of rj on the behaviour of bouncing solutions 
in Fig. m for both closed and open universe. All the plots are generated for the same set of 
initial conditions and the same set of parameters Wm = 1/3, sign{pk) = +ve but with three 
different values of parameters r] =1/4, 1/6 and 1/8. It has been observed that the value of p 
has a direct impact on the occurence of bouncing point. Indeed, the position of bouncing point 
is delayed as we decrease the value of p for both closed and open universe as shown in Fig. [TO] 
{top left and top right). The bottom left and bottom right of Fig. [TO] indicate the effect of p 
on the curvature parameter for both closed and open cases respectively. It has been found that 
the magnitude of maximum value of z, at the bounce, decreases as we decreases the value of p. 
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Figure 7: tt = ^ vs x for closed {left) and for open {right) with Wk = —2.0, Wm = 1/3 and 
sign{pk) = +ve 




Figure 8: z vs x for closed {left) and for open {right) with Wk = —2.0, Wm = 1/3 and sign{pk) = 
-\-ve 


X 




Figure 9: Comparision between closed and open for x vs N {left) and Qm vs N {right) with 
Wk = —2.0, Wm = 1/3 and sign{pk) = +ve 
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Figure 10: Evolution x for closed {top left) and open {top right); z for closed {bottom left) and 
open {bottom right) with values rj = 1/4,1/6, and 1/8, Wm = 1/3 and sign{pk) = +ve 


13 



















5 Conclusion 


A cosmological scenario with a noncanonical scalar field and matter is explored in this work. 
Using dynamical equations for a set of dimensionless dynamical variables, we find all fixed point 
for the two cases with positive and negative kinetic energy density term in FRW closed and open 
universe. Allowed region of parameter spaces for stability of fixed points are shown for both 
cases. The necessary and sufficient conditions for a nonsingular bounce are obtained in terms 
of the dynamical variables. Thus, stable bouncing solutions are obtained satisfying nonsingular 
bouncing conditions and stability criteria. This is achieved for the negative energy density of 
matter, Clm, with equation of state parameter Wm = 1/3 in both closed and open universe. In 
addition to this, the finitude of curvature parameter at the bounce is obtained as expected in 
a nonsingular bouncing scenario and universe becomes flat at late time irrespective of whether 
we start with a closed or open one. Finally, the effect of the parameter rj on the behaviour 
of bouncing solution is noted. It is seen that the point of occurence of bounce is delayed as 
we decrease the value of rj and the magnitude of value of curvature parameter at the bounce 
decreases with ij for both open and closed universe. 

We restrict our analysis to a positive sign of potential. It is straightforward to extend our 
analysis for a negative potential by changing the parameter sign{y) to —1. 
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